Wind reversals in turbulent Rayleigh-Benard convection 
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The phenomenon of irregular cessation and subsequent reversal of the large-scale circulation in 
turbulent Rayleigh-Benard convection is theoretically analysed. The force and thermal balance 
on a single plume detached from the thermal boundary layer yields a set of coupled nonlinear 
equations, whose dynamics is related to the Lorenz equations. For Prandtl and Rayleigh numbers 
in the range 10 -2 < Pr < 10 3 and 10 7 < Ra < 10 12 , the model has the following features: (i) 
chaotic reversals may be exhibited at Ra > 10 7 ; (ii) the Reynolds number based on the root mean 
square velocity scales as Re rm3 ~ Ra' ' 41 ' °' 47 ' (depending on Pr), and as Re rms ~ p r -[o.66 ■■• 0.76] 
(depending on Ra); and (iii) the mean reversal frequency follows an effective scaling law u/(y L~ 2 ) ~ 
p r -(o.64±o.oi) p^a - 44 * - 01 The phase diagram of the model is sketched, and the observed transitions 
are discussed. 



One important issue in turbulent Rayleigh-Benard 
convection is the interplay between the large-scale circu- 
lation (the so-called wind) 0] and the dynamics of plumes 
detached from the thermal boundary layers Q. In par- 
ticular, such interplay seems to be relevant in the process 
of circulation reversals .which occur in an irregular time 
sequence 0, 0, IE IS Hi ■ Remarkably, similar reversals 
are also observed in the wind direction of the atmosphere 
and in the magnetic polarity of the earth ^(j . 

In principle, two reversal scenarios are possible: Rever- 
sal through cessation of the convection roll, and reversal 
through its azimuthal rotation. With two temperature 
sensors placed close to each other near the sidewall 0, |^ , 
one can detect roll reversals, but not distinguish between 
the two scenarios. With several sensors placed along the 
azimuth of the cell, Cioni et al. Q succeeded to detect 
reversal through azimuthal rotation of the roll. Reversal 
through rotation was also detected in refs. 0,0. How- 
ever, with an ingenious multi-probe setup, Brown, Niko- 
laenko, and Ahlers Q were able to distinguish between 
the rotation and cessation scenarios, and many reversals 
through cessation were detected. Reversal through ces- 
sation was also observed in two-dimensional numerical 
simulations of the Boussinesq equations (see fig. 8 of ref. 
[lll | and fig. 12 of ref. 0])> where the rotation scenario 
is of course impossible. 

Since reversal through cessation is a more surprising 
scenario, the aim of the present work is to reveal its phys- 
ical mechanism. Qualitatively, the picture is as follows 
[l3j : If an uprising hot plume gets too fast because of 
a temperature surplus, it fails to cool down sufficiently 
when passing the top plate. It then is still warmer than 
the ambient fluid when advected down along the side- 
wall. By buoyancy it therefore looses speed and coun- 
teracts the large-scale circulation. Indeed, the downward 
wind may be counteracted so strongly that it stops or 
even reverses its direction. This mechanism can be effec- 
tive only for sufficiently strong wind, i.e., for sufficiently 



large Reynolds number, because for slow motion the ther- 
mal diffusivity n has enough time to reduce the tempera- 
ture surplus of the originally warmer plume relative to its 
neighbourhood. Then its power to reverse the circulation 
by buoyancy is gone. 

The model: In order to quantify the cessation mecha- 
nism discussed above, let us first characterize the size of 
a circulating plume. As shown in figure Q a single plume 
will be understood as a thermal structure of width Ag 
(the thickness of the thermal boudary-layer from which 
it originated) and length L (the height of the convection 
container). In addition, its volume is assumed to scale 
as Xg L, with a typical cross-section area A 2 ,, and surface 
area A« L. 
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FIG. 1: Sketch of the motion of a single plume of width As 
and height L. In an aspect-ratio-one container, the circulation 
radius is given by L/2. 



Supposing that such a plume circulates with velocity 
U(t), it is reasonable to expect that its dynamics is essen- 
tially a matter of balance between buoyancy and drag. 

In the Boussinesq approximation the buoyancy force 
(per mass) is given by = — a p (T — T )g, where 
a p is the isobaric thermal expansion coefficient, T the 
plume temperature, To the mean temperature, and g the 
gravitational acceleration. On the other hand, the drag 
force (per mass) on the plume has the strength fd — 
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\ C(Re) U 2 L~ x , where C(Re) is the drag coefficient, and 
the Reynolds number is defined by Re = LU 1/ . Here, 
C(Re) is taken as |T^ : 



C(Re) 





8Re^ 



(1) 



where b = 8.4 is the Kolmogorov constant. Equation 
||TJ describes the transition from the strongly decreasing 
drag ~ 1 /Re in the viscous regime to the Re- independent 
drag in the turbulent regime. As pointed out in ref. |l4| . 
it pretty well agrees with experimental data. 

Now, let us consider the thermal interaction between a 
single plume and its surrounding. Strictly speaking, the 
surrounding consists of the fluid as well as the sidewalls, 
and the top and bottom plates. We do not distinguish 
between all these and describe the temperature of the 
plume surrounding T s (ip) by a time- independent profile: 



T s (<p) = T + y cos if, 



(2) 



where A is the temperature difference between the hori- 
zontal plates. We Fourier expand the temperature vari- 
able of the plume: 

oo 

T(<p,t) = T + Y, [An(t) sin(n^) + B n (t) cos(r^)] , (3) 

n=l 

where A n (t) and B n (t) are the amplitudes. 

Equations of motion: In order to derive the equations 
of motion for a single plume, we follow an analogy with 
the Malkus waterwheel 0, rLfij . On the basis of this 
analogy, our intent is to acquire an understanding of the 
wind dynamics through nonlinear model equations. 

To begin, let us consider the balance of forces (per 
mass) on the plume: 



dU_ 
~dt 



fb(<P,t) simp - f d . 



(4) 



Substituting the previous relations into Q , and integrat- 
ing the resultant expression with respect to tp from to 
2-7T, one readily finds: 



dU 1 
dt 2 py 



1 U 2 



(5) 



Remarkably, the temporal behavior of U is coupled to 
the amplitude of the first temperature mode A\ only. 

The temporal change of the plume temperature is given 
by advection and by diffusion. For the latter, we assume 
a relaxation ansatz for the temperature deviation T — T s 
from the surrounding, with the diffusive time scale t k = 
XqL/k, i.e., 



dT 
~dt 



U dT 
1/2 £hp~ 



T-T K 



(6) 



The physics behind the definition of r K is that the thermal 
loss is proportional to the plume surface, and inversely 
proportional to the thermal diffusivity. 

Substituting (0 and J3J into ©, and equating the 
coefficients of each harmonic separately, one obtains: 



dAi 
~dT 
dB x 
~dT 



K B + kA 2 jja (8) 

XgL 2\qL L 



We write the three coupled ODEs JSJl, J7J), and ||SJ| in 

nondimensional form. The dimensionless variables are 
X = 2~ 1 Nu" 2 k" 1 LU, Y = A x , Z = (Ra" 1 - 

2_BiA _1 )r, r = 2NukL -2 £, and the dimensionless con- 
trol parameters read: 



9 Pr 
4 Nu 



1 Ra 



18Nu Ra c 



(9) 



where Ra c = 1708, Pr = v/k is the Prandtl number, and 
Ra = a v gl? A I [y k) the Rayleigh number. The Nusselt 
number Nu comes from the relation Xg/L = 1/(2 Nu). 
Then, the system of equations J7J and © becomes: 



dX_ 
~dV 

dY 

~aW 
dZ_ 



oY -aX 



= rX - Y - X Z, 



= -Z + XY. 
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(10) 

(11) 
(12) 



This system resembles the Lorenz equations 
which have also been used to describe laminar flow con- 
fined in a toroidal loop 0,H(J. Here equations JTDJl-lO 
have been derived to model plume reversals in the tur- 
bulent regime. They will be referred to as the modified 
Lorenz equations. There are two essential differences as 
compared to the standard Lorenz system |2jj: (i) The 
parameters a and r are related to the Nusselt number, 
which is known to follow a nonuniversal (Pr-dependent) 
scaling with Re j2^|. This is a key difference, since in 
the Lorenz equations a = Pr and r = Ra/Ra c . (ii) The 
ordinary differential equation for X has a new nonlinear 
term, due to the turbulent drag on the plume. 

Phase diagram: To investigate the dynamical proper- 
ties of the system {TUJ-lO, we have scanned the pa- 
rameter space Ra x Pr in the range 10 7 < Ra < 10 12 , 
10~ 2 < Pr < 10 3 . Technically, our numerical scheme 
was based on a fourth-order Runge-Kutta method [22| . 
with adaptive stepsize control in time, and increments 
of 0.1 for log 10 (Pr) and log 10 (Ra). The Nu-input re- 
quired for coefficients 10 is provided by Grossmann- 
Lohse theory [2^, and as initial condition we adopted 
(X = 1, Y = 1, Z = 1). As for check with other initial 
values see later. 
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FIG. 2: Time series of the dimensionless plume velocity X, 
for Pr = 0.1 at (a) Ra = 10 10 ' 8 (uniform circulation), (b) 
Ra = 10 10 ' 9 (chaotic reversals), and (c) Ra = 10 110 (periodic 
reversals) . 
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FIG. 3: Phase diagram in the Ra x Pr plane: at sufficiently 
large Ra, the state of uniform circulation (blank region) gives 
place to chaotic (stripped region) or periodic (dotted region) 
wind reversals. Note the small periodic windows in the chaotic 
range. 

An insight into the structure of the phase diagram can 
be acquired by considering some representative times se- 
ries of X(t). In particular, for fixed Pr and increasing Ra, 
three examples are shown in figure|2l first, a state of uni- 
form circulation [cf. plate (a)] ; then emergence of chaotic 
reversals [plate (b)]; and, ultimately, periodic reversals 
[plate (c)]. Figure shows the phase diagram in Ra x Pr 
space, displaying a sharp onset between the steady and 
the reversal domain. We emphasize that the transition 
curve between these domains remains unchanged for a 
variety of initial conditions. 

Onset of reversals: The onset of reversals can be un- 
derstood in terms of the typical time scales of the sys- 
tem: the thermal diffusion time r K and the turnover 
time tjj = ttL/(U), where (•) denotes the time aver- 
age. Qualitatively, it is reasonable to expect wind re- 
versals when tjj <C t k , because in such case the cir- 
culation is so fast that the plume has no time to lose 
its temperature contrast. Indeed, we find that the ratio 
t u/ t k — 27rNuPr~ 1 (Re) -1 is a monotonically decreasing 
function of Ra for constant Pr, roughly proportional to 
Ra -1 / 6 . The overall form of the onset curve well resem- 
bles its counterpart in the phase diagram of the Lorenz 
model cf. Dullin et al. El. 



Reynolds number: We now come to the dependence 
of the variance of the Reynolds number Re rms = 
LurmsV^ 1 based on the root mean square velocity 
Ums = \J ((U - (U)) 2 ). Figure H shows Re rms (Ra, Pr): 
In plate (a), the Ra-scaling exponent increases from 0.41 
to 0.47 for increasing Pr from 0.7 to 316; in plate (b), 
the Pr-scaling exponent decreases from —0.66 to —0.76 
for falling Ra from 10 12 to 10 9 . Experimentally, a similar 
Pr-dependence has been reported |2^| for the Reynolds 
numbers based on the maximum wind velocity, on the 
oscillation frequency of the large-scale circulation, and 
on the rms velocity. 
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FIG. 4: Reynolds number variance Re rms based on the root 
mean square velocity as function of (a) Ra and (b) Pr. 

Mean reversal frequency: The abrupt change of X(t) 
with r [cf. figure Efb)] suggests that the wind switching 
can be approximately considered as an almost instanta- 
neous event represented by the moment at which it oc- 
curs. Here, we follow Sreenivasan et al. Q and define t' n 
as the interval between an arbitrary origin in time and 
the nth wind reversal. Similarly as in |4|, we also find a 
linear relation t' n ~ n, which suggests a mean interval t' 
between reversals. In this way, we define u> = 1/t' as the 
mean reversal frequency, and its dimensionless counter- 
part as Q = u L 2 /v. 
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FIG. 5: Mean reversal frequency as function of (a) Ra at 
fixed Pr and (b) Pr for given Ra. 

FigureElshows that u ~ Pr -(°-65±0.oi) Ra o.44±o.oi Tq 
our knowledge, the only experimental measurement of 
the reversal frequency has been carried out in cryogenic 
helium gas [ 



namely Cj ~ Ra for Pr = 0.75 and 
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2.1 x 10 8 < Ra < 1 x 10 13 . The disagreement be- 
tween our result and the particular measurement sug- 
gests that a model based on only 3-modes for the plumes 
is quantitavely inadequate. In this qualitative sense, our 
simple deterministic system well mimics the dynamics 
of reversals, and is a complementary approach to the 
stochastic model of noise- induced switchings between two 
metastable states 0. Here, the Lorcnz attractor itself 
captures the bistable transitions, but a more quantita- 
tive description of the reversal phenomenon (also includ- 
ing the rotation scenario) would involve a subtle combi- 
nation of deterministic chaos and noise. This could be 
done in the spirit of ref. [13] (sec. Ill C), and is left for 
future work. 

We thank G. Ahlers, K. R. Sreenivasan, and J. Niemela 
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